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Abstract. 

The contributions to the g — 2 of the muon from some eighth-order (four- loop) graphs 
containing one-loop and two-loop vacuum polarization insertions have been evaluated 
analytically in QED perturbation theory, expanding the results in the ratio of the electron 
to muon mass (m e /m M ). The results agree with the numerical evaluations and the 
asymptotic analytical results already existing in the literature. 
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Recently, the sixth-order (three-loop) coefficient of the QED perturbative 

expansion of the difference between muon and electron (g-2) 
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was calculated in closed analytical form [1][2]. The eighth-order (four-loop) coefficient 
affi — is known only in numerical form [3] ; its more recent value is [4] : 



ajf) - c4 4) = 127.55(41) . 
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In this work we have calculated in analytical form the contributions to ajL — ai^ from 
some graphs containing insertions on a single photon line of one-loop and two-loop vac- 
uum polarization subdiagrams (see fig.l); the considered graphs are shown in figs. (2) and 
(3). 

We found convenient to express the results expanding them in the ratio of the 
electron and muon masses (m e /m M ); for the sake of extensive numerical checks we have 
calculated the terms containing up to (m e /m M ) 16 . Unfortunately the coefficient of each 
term becomes more and more cumbersome as the power of (m e /m^) increases, so that 
we will list here the terms of the expansions up to (m e /m M ) 2 only. 

The analytical expressions of the contributions to the muon anomaly of the graphs 
shown in figs. (2) and (3), accounting for the proper multiplicity factors, are (r = m e /m M ): 
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Here ((p) is the Riemann ^-function of argument p, ((p) = — , (whose first values 

n=l 

are C(2) = tt 2 /6, C(3) = 1.202056903..., C(4) = tt 4 /90, C(5) = 1.036927755...), a 4 = 

oo ^ oo ^ 

V 7 = 0.517479061... , a 5 = V F = 0.508400579... , and (3 2 is the Catalan 

^— ' 2 n n 4 2 n n° 

n=l n=l 

constant (3 2 = V 777 = 0.915965594... 

H ^ 2n + l 2 

We note the appearance of transcendental constants (3 2 in eq.(6) and a$ in eq.(8) and 
eq.(9); these constants appear for the first time in a contribution to the lepton anomaly. 

Of the above expressions, eqs.(3)-(ll), only the leading terms of eq.(3) and eq.(6) 
were already worked out in analytical form using renormalization group techniques [3] [5] . 
( 1 ). An examination of eqs.(3)-(ll) shows that the expressions of the contributions of the 
graphs containing only electron loops, eqs.(3),(6),(9) and (11), have a linear {m e /m^) 
term; the numerical values of the corresponding coefficients are respectively 6.41, -5.61, 
-0.44, -2.84. The (m e /m^) expansions of the contributions of the other graphs (containing 
also muon loops) begin with the (m e /m M ) 2 term. 

We list now the numerical values of eqs.(3)-(ll) obtained using the experimental 
value [7] (m^/m e ) = 206.768262(30) and taking into account all the calculated terms of 
the expansions in {m e /m^), compared with the corresponding numerical values given in 



( 1 ) We want point out that, contrarily to the assertion of Appendix A of ref.[3], even 
the leading terms of the sum of the contributions from the graphs of fig. 2(e) and fig. 3(c) 
can be calculated in analytical form using renormalization group techniques. In fact, 
the integral I\ of eq.(A.25) of ref.[3], evaluated numerically in that work, can be worked 
out in analytical form using the contributions of sixth-order graphs containing vacuum 
polarization insertions which are known in analytical form by long time (see eq.(2.22) 
of ref. [6]); the analytical value of this integral is found to be Ii = —it — -tt In 2 — 

^ In 4 2 - 8a 4 - 8C(3) + ^tt 2 In 2 - 4tt 2 + ^ which, inserted in eq. (A.26) of ref. [3] , gives 
the leading terms of the sum of our eq.(7) and eq.(ll). 
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ref.[3]: 
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The total sum is 

ajf) [fig.2 + 3] = 3.2048378(8) 



(our); 



3.1845(66) (ref. [3]). 
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The numerical error of our results is induced by the experimental uncertainty of 
(m M /m e ); in order to reach such a precision accounting of terms up to (m e /m M ) 4 is 
needed. Note that the linear terms in (m e /m M ) are essential to check the results of ref. [3] 
within their precision: as an example, the contribution of the linear term of eq.(3) is 
0.031 which is about 23 times the error of the correspondent numerical result of ref. [3]. 
We found that our results are in good agreement with numerical results of ref. [3]; only 
the contribution affi [fig.3(a + b)] shows a slight disagreement, at the level of 3.6cr. This 
is a remarkable cross check of results, due to the difference of the methods followed in 
the two derivations. 

Let us now consider the contributions to the electron anomaly from the graphs shown 
in figs. (2) and (3) with \i and e leptons exchanged ( 2 ); we list only the leading terms of 
the expansions (r = m e /m M ): 
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( 2 ) The contributions to the electron anomaly of the graphs shown in figs. (2) and (3) 
with \i leptons replaced by electrons can be found in ref. [8] . 
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due to the smallness of the ratio (m e /m M ), these contributions are almost negligible, the 
numerical value of the sum of eqs.(13)-(21) being 

c4 4) [fig.2 + 3] = -1.796 x 10" 4 . (22) 



value for a^P — ai^ 



Using eq.(2) and our results (12) and (22) we can work out a new slightly different 

ajf) - a { e 4) = 127.57(41) . (23) 
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We sketch the method used for obtaining the analytical expressions of eqs.(3)-(ll) 
and eqs.(13)-(21). Quite in general the contribution to the anomaly of the muon from a 
vertex graph with vacuum polarization insertions with electron loops can be written as 




(24) 




where ImII n I — ^ ) * s the imaginary part of the vacuum polarization insertion in nth 



order and K\ I — - is the anomaly contribution from some set of /th-order vertex graphs 



in which a photon line has been given a mass b. Both quantities are analytically known 
up to fourth order [6] [9] [10]. 

Once that the suitable expressions of K and Imll are inserted in eq.(24), the contri- 
bution to the muon (g-2) becomes a sum of one-dimensional integrals containing square 
roots, logarithms, dilogarithms and at worst trilogarithms of the variable b. 

Unlike ref.[2], where analogous sixth-order integrals were calculated in closed an- 
alytical form, we found convenient to calculate the integrals of eq.(24) expanding in 
(m e /m M ). We split the integration region into two parts by introducing a cut A, such 
that ml C A < m\ 2 : in the region where b < A we find C 1 so that we can expand 



Ki ( — - ) for small values of the argument, while in the region where b > A we expand 



Imil n — - for large — ~ . Integrals over these regions were calculated analytically 



using the method described in ref.[2] as functions of m e , m M and A; summing up the 
analytical contributions of the two regions the dependence on A drops out, as expected. 
A similar method is used when the e and \x leptons are exchanged. 

Finally, we checked that the direct numerical evaluations of the integrals (24) are in 
perfect agreement with our analytical expressions for the same quantities. 
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Figure captions 



Fig.l: The irreducible fourth-order vacuum polarization subdiagram. 

Fig. 2: Eighth-order vertex graphs obtained with insertions of second- and fourth-order 
vacuum polarization subdiagrams on the second-order vertex graph. 

Fig. 3: Examples of eighth-order vertex graphs obtained with the insertion of second- and 
fourth-order vacuum polarization subdiagrams on fourth-order vertex graphs. 
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